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Abstract. New optical photometric and spectroscopic data of the bright hypernova SN 2002ap are presented. 
The obtained BVRI light curves as well as the optical spectra agree well with other published data series. The 
distance has been inferred by applying the Expanding Photosphere Method for the data around maximum. The 
derived 6.7 Mpc is in good agreement with recent photometric distances of M74. However, the total (random plus 
systematic) uncertainty of the EPM-distance is at least ±4.5 Mpc (about 70 %). The physical parameters of the 
SN have been determined via simplified analytic models of the light and velocity curves. It is confirmed that SN 
2002ap was a less energetic hypernova, the kinetic energy was 4 - 8 ■ 10^^ erg, and the reddening-free absolute 
bolometric magnitude reached —16.63 mag (corresponding to Lboi ~ 3.47 ■ 10*1/0 ), about 2 magnitude less than 
the prototype hypernova SN 1998bw. 
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1. Introduction 

■ The idea that core-collapse supernovae (SNe) may pro- 
duce observable gamma-ray bursts, has received many the- 
oretical as well as observational attention since the discov- 
ery of SN 1998bw w ithin the error box of GRB 980425 
l|Galama et al. 1998|l . Recently, the optical afterglow of 
GRB 03032 9 that showed SN 1998bw-Uke spectrum 
(SN 2003dh,l(:^arnavich et al. 20031 IT^hornock et ai. 20031 
IStanek et al . 20031 |Lipkin et al., 2004| ) provided a very 
convincing evidence of the supernova - GRB connection. 
With GRB 031203 - SN 20031w IjGal-Yam et al. 2004|l 
there are now three low-redshift (z < 0.3) GRBs that 
are all associated with SNe. The spectra of these SNe 
were characterized by shallow, very broad lines and a 
very blue continuum at the early phases. Usually, their 
spectroscopic classification is Ic, although they do not re- 
ally resemble other Type Ic SNe. The very broad lines 
in the early spectra are indicative of high expansion ve- 
locity, up to 40,000 kms"^ Such SNe are now referred 
as "hypernovae" , because of their high expansion veloc- 
ity as well as large kinetic energy released by the ex- 
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plosion (llwamoto et al. 1998|l . Actually, some model com- 
putations ( Hoflich, Wheeler fc Wang, 1999| ) showed that 
such high expansion velocities can also be consequences of 
an asymmetric explosion seen at different viewing angles 
rather than extremely high intrinsic kinetic energy. 

SN 2002ap, discovered by Y. Hirose on Jan. 29, 
2002 (|Nakano et al. 2002|l . turned out to be a very 
interesting event in many respects. First, it occured 
in M74 (NGC 628), a nearby spiral galaxy (d « 7 
Mpc), thus, it is one of the closest SNe observed ever. 
Second, it was discovered well before maximum, and 
the early phase spectra were very similar to those of 
SN 1998bw. Th erefore, SN 2002ap was quickly classified 
as a hypernova l|Meikle et al. 20021 |Kinugasa et al. 2002| 
|Gal-Yam fc Shemmer, "2002) 1. Third, despite the expec- 
tations, SN 2002ap could not be associated with 
any detected GRBs UTotani, 2003| ). Some authors 
l|Leonard et al. 2002a|l suggested that SN 2002ap should 
be classified as "Peculiar Type Ic" SN instead of hyper- 
nova, because in this case the latter term can be reserved 
to those SNe that actually produce GRBs. 

In spite of the lack of associated GRB, 
SN 2002ap has been extensively foUowed-up in many 
bands extending from radio ( |Berger et al., 2002) ) 
to X-ray IjSutaria et al. 20021 ISoria et al. 2003|l . 
Optical spectroscopy and/or photometry were 
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reported in many papers fGal-Yam et al. 2002" 
Borisov et al. 20021 [Cbok et al. 2002. .Pandey et al. 2002 
Kinugasa et al. 2002| lYoshii et al. 20031 

Fble^'eTlir'^ OOS^ Based on optical spectrophotom- 
etry, significant continuum polarization was de- 
tected (jKawabata et al. 20021 ILeonard et al. 2002al 
|Wang et al. 2002| ), suggesting an asymmetric explosion. 

From these observations, a more-or-less consistent 
physical picture was revealed, containing a relatively 
low-mass progenitor (M = 2 — 5 Mq), exploding with 
E ^ A - 10 10^^ ergs kinetic energy, synthesizing 0.07 
- 0.1 Mq of ^^Ni Ipalzah et al. 20021 IMaeda et al. 20031 
[Foley et al. 2003] ). The original mass of the progenitor was 
estimated as Mzams ~ 20 — 25 AIq. However, the lack 
of the detection of the progenitor on previous optical and 
near-IR frames of M74 fSma rtt et al. 2002jl may imply a 
lower mass and/or lower luminosity. 

All of these physical parameters are based on the 
assumed distance to the host galaxy of SN 2002ap. 
There are two modern distance estimates of M74 in 
the literature. ISharina et al. 19961 derived /io = 29.32 
mag as a reddening-free distance modulus, based on 
the luminosity of the brightest blue supergiants. The 
Sharina et al. result confirms the distance determined 
by 'Sohn fc Davidge, 1996| (fip = 29.3), although Sohn & 
Davidgc claimed that blue supergiants are unreliable dis- 
tance indicators and they preferred the use of the brightest 
red supergiants for this purpose. However, both groups as- 
sumed the total absorption in i? as As = 0.13 mag based 
on the map of [Burstein fc Heiles, 19821 The more recent 
reddening map oflSchle gel et al. 1998| predicts = 0.301 
mag in the direction of M74, which reduces the above dis- 
tance moduli to /io = 29.15 mag, corresponding to d = 6.8 
Mpc. 

Surprisingly, the various attempts to determine the 
distance of M74 have led to very inconsistent re- 
sults. For example, [Sandage fc Tammann, 1974| derived 
(Uo = 31.46 mag from the sizes of Hll-regions, while 
IBottinelh et al. 19841 inferred = 26.60 from Tully- 
Fisher relation. Thus, there is an order of magnitude dis- 
agreement between the distances determined by different 
methods. 

In all papers cited above, the distance of SN 2002ap 
was assumed to be 7.3 Mpc, corresponding to the original 
distance modulus of ISharina et al. 19961 This value itself 
may introduce a 0.2 mag systematic uncertainty in the 
derived absolute magnitudes of the SN (which then incor- 
porates into the other inferred physical parameters), be- 
cause this distance is based on the earlier reddening map. 
Moreover, the large scatter in the early M74 distances 
may further complicate the situation. A new, independent 
distance measurement of M74 based on SN 2002ap itself 
seems to be useful to clarify the picture. 

In this paper we present a new distance estimate of SN 
2002ap and M74 via the Expanding Photosphere Method 
(EPM). For this purpose we obtained new optical spectro- 
copic and photometric measurements of SN 2002ap, from 
t = —A days up to t = -f291 days. These are described in 
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Fig. 1. SN 2002ap and the local comparison stars in the 
field of M74. 



Section 2. The distance estimate is discussed in Section 3. 
Section 4 contains the analysis of the bolometric light and 
radial velocity curves, based on a simple analytic model. 
Section 5 summarizes our results. 

2. Observations and data reduction 

2.1. Photometry 

The photometric observations were started on Jan 
31, 2002, just after the discovery of SN 2002ap. 
Most of the data were obtained with the 60/90 cm 
Schmidt-telescope of Konkoly Observatory (installed 
at Piszkesteto Mountain, Hungary) equipped with a 
1536x1024 Photometries camera. Some CCD frames were 
also collected using the 1 m RCC telescope at Konkoly 
and the 40 cm Cassegrain at Szeged Observatory (see 
IVinko et al. 20031 for the technical details of these instru- 
ments) . 

The CCD frames were reduced using the standard way 
in IRAF^. Comparison stars were selected in the vicinity 
of the SN (Fig.l), and differential magnitudes were deter- 
mined via aperture photometry. Because the SN occured 
outside the visible region of its host, the background was 
not affected significantly by the light of the host galaxy. 
Thus, aperture photometry could be easily applied. 

The instrumental magnitudes were transformed into 
the standard system using the equations given by 
IVinko et al. 20011 The comparison stars were calibrated 
via standard stars in the field of M67. Then, the standard 
magnitudes were compared with those given by Henden 
l|Oal-Yam et al. 2002|l and Foley et al. (2003). The three 
datasets showed acceptable agreement within 0.03 mag 

^ IRAF is distributed by NOAO which is operated by the 
Association of Universities for Research in Astronomy (AURA) 
Inc. under cooperative agreement with the National Science 
Foundation 



J. Vinko et al: The distance of SN 2002ap 



3 



Table 1. Final adopted magnitudes of local standard 
stars. Errors are given in parentheses. 



Star 


V 


B-V 


V -R 


V -I 


A 


12.91(0.05) 


1.24(0.06) 


0.78(0.02) 


1.46(0.04) 


B 


13.80(0.02) 


0.72(0.01) 


0.42(0.02) 


0.81(0.03) 


C 


13.08(0.03) 


0.76(0.02) 


0.45(0.01) 


0.89(0.03) 


D 


13.24(0.02) 


0.77(0.02) 


0.44(0.01) 


0.87(0.02) 




Fig. 2. Comparison of the standard sequences of Henden 
and Foley et al. plotted as a function of our standard mag- 
nitudes. Left panel: V magnitudes (open circles: Henden, 
filled circles: Foley et al.), right panel: B — V (circles), 
V — R (squares) ,V — I (triangles) . 



V) = 0.09 mag (see next section), the magnitudes were 
converted into monochromatic fluxes using the calibration 
of lBesseU, 2001) The near-IR contribution was estimated 
from the JHK data of lYoshii et al. 2003L while contribu- 
tion from other bands were neglected. Then, the bolomet- 
ric fluxes were calculated by integrating the monochro- 
matic fluxes with a simple trapezoidal rule. This resulted 
in the "observed" bolometric flux as a function of time. 
Finally, the absolute bolometric magnitudes have been de- 
termined applying the true distance modulus fiQ = 29.13 
mag derived in this paper (see Sect. 4). 

The absolute bolometric magnitudes have been com- 
pared with those derived by lYoshii et al. 20031 Because 
Yoshii et al. applied /zq = 29.5 mag as distance modu- 
lus, their data were corrected to match the slightly lower 
distance modulus used above. After correcting for this dif- 
ference, the two bolometric light curves showed very good 
agreement. Thus, it is concluded that the various deter- 
minations of the UVOIR bolometric light curve of SN 
2002ap have led to consistent results, and the source of 
the main uncertainty of the bolometric light curve is the 
true distance modulus of the SN. 

The bolometric light curve of SNe can be described 
analytically by the simple model of IContardo et al. 20001 
In this model the light curve has the form 



mo + j{t - to) + go cxp[-{t - tp) V(2crg)] 



(Fig. 2), except for the brightest star (A), where the differ- 
ences were slightly higher. The final adopted magnitudes 
of our comparison stars are determined as the average 
of the three standard sequences. These are collected in 
Table 1. The uncertainties are estimated from the stan- 
dard deviation of the averaged data. 

The magnitudes listed in Table 1 were used to cal- 
culate the calibrated standard magnitudes of SN 2002ap. 
The final results (after averaging the magnitudes obtained 
from star A, B, C and D) are listed in Table 2. Fig. 3 
shows the light curve in BVRI bands. In order to com- 
pare our photometry with other published data, we also 
plot photometric data of SN 2002ap collected from liter- 
ature. The agreement is satisfactory (note that the tabu- 
lated light curve data in Foley et al. (2003) already con- 
tain the reddening correction - this was removed before 
plotting in Fig. 3). As usual, the phase is referenced to 
the moment of i3-maximum, which was determined as JD 
2452311.5 (note that it is 2 days less than the one given 
hvl(;aJ-Ya,m et al. 200211. 

2.2. Constructing the bolometric ligtit curve 

The bolometric light curve was constructed after combin- 
ing all available optical and near-IR measurements col- 
lected from the literature. After dereddening with E{B — 



Note that we applied only one gaussian component, 
because the bolometric light curve of SN 2002ap does not 
show significant bump on the descending branch. The light 
curve parameters were determined via least-squares fitting 
as mo = 13.82 ± 0.1, to = 310.0 ± 0.5, t2 = 297.95 ± 
0.02, 7 = 0.0172 ± 0.0005, ctq = 18.3658 ± 0.0004, era = 
3.1677 ±0.0002, go = 1.4684 ±0.2. Among these, the most 
interesting parameter is the late-time decline rate 7 that 
is physically linked to the Co-Fe decay and the radiation 
transport in the expanding SN ejecta. The value derived 
above, 7 = 0.017 mag/day is in good agreement with the 
one published by |Pandey et al. 2002| (0.0199 ± 0.0004). 

In Fig. 4 we plot the "observed" bolometric light curve 
(wboi = Miioi + ^0) together with the analytic model 
given above and the expected decline rates due to the Ni- 
Co-Fe radioactive decay. As seen in Fig. 4, the decline of 
SN 2002ap is significantly faster than the Co-Fe radioac- 
tive decline rate (r = 0.0098 mag/day), as also noted by 
IPandey et al. 2002| and lYoshn et al. 20031 They explained 
this as a leakage of 7— photons from the atmosphere, indi- 
cating a transparent, less massive ejecta. Our new late- 
time photometry fully supports these conclusions. The 
connection between the physical state of the expanding 
atmosphere and the bolometric light curve is analyzed fur- 
ther in Section 4. 
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Fig. 3. The liglit variation of SN 2002ap around maximum (left panel) and during the nebular phase (right panel). 
The data obtained with different filters were shifted vertically for better visibility (with -f 0.5, 0.0, —0.5 and —1.0 mag 
for B, V, R and / data, respectively). The symbols denote the following data sources: plus signs: Foley et al. (2003); 
crosses: Pandey et al. (2002); open squares: Yoshii et al. (2003); filled squares: Gal- Yam et al. (2002); asterisks: present 
paper. 



Table 2. Photometry of SN 2002ap. The telescope codes are: A - 60 cm Schmidt, Konkoly, B - 1 m RCC, Konkoly, 
C - 40 cm Cassegrain, Szeged. 



JD 


B 


V 


R 


I 


Tel. 


52307.30 


13.27(0.013) 


12.79(0.008) 


12.89(0.009) 


13.02(0.008) 


A 


52308.30 


13.13(0.015) 


12.62(0.006) 


12.71(0.008) 


12.87(0.007) 


A 


52316.30 


13.35(0.080) 


12.53(0.043) 


12.34(0.031) 


12.40(0.032) 


C 


52322.30 




12.91(0.041) 


12.55(0.022) 


12.39(0.024) 


C 


52323.30 




12.89(0.067) 




12.45(0.035) 


C 


52325.30 




12.99(0.029) 


12.74(0.049) 


12.53(0.042) 


C 


52328.30 


14.43(0.051) 


13.33(0.046) 


12.87(0.052) 


12.59(0.028) 


C 


52337.20 


15.05(0.010) 


13.99(0.007) 


13.48(0.008) 


13.01(0.009) 


A 


52338.20 




14.01(0.014) 


13.51(0.008) 


13.07(0.015) 


A 


52339.20 


14.96(0.029) 


14.10(0.016) 


13.59(0.021) 


13.08(0.019) 


A 


52481.40 




17.18(0.007) 


16.19(0.008) 


16.14(0.006) 


A 


52482.40 




17.16(0.006) 


16.25(0.005) 


16.08(0.013) 


A 


52484.40 




17.03(0.015) 


16.28(0.033) 


16.08(0.050) 


A 


52526.50 


18.32(0.022) 


17.80(0.015) 


16.78(0.008) 


16.83(0.012) 


A 


52514.60 




17.61(0.027) 


16.55(0.030) 


16.65(0.018) 


B 


52602.40 




18.94(0.190) 


18.35(0.110) 


17.80(0.125) 


A 



2.3. Optical spectroscopy 

The spectroscopic observations were carried out at David 
Dunlap Observatory (DDO) between the 2nd and 25th 
February, 2002 {t = —3 to +20 days), with the Cassegrain 
spectrograph attached to the 74" telescope. Two different 
setups have been applied: the 100 lines/mm grating was 
used for the low dispersion observations, while the 1800 
lines/mm grating was selected for obtaining high resolu- 
tion spectra. The reciprocal dispersion was 3.60 A/pixel 
and 0.2 A/pixel, respectively. The low-resolution spectra 
extended from 4000 to 7500 A, while the high-resolution 
spectrum was centered on Na D and covered a 200 A- wide 
spectral interval. The slit width was 306 microns (cor- 
responding to 1.8 arcsec on the sky) for both observing 



modes. The seeing was between 2-3 arcsec during the 
observations. 

The reduction was done in the standard way using 
IRAF. Flux calibration was performed via the spectra of 
comparison stars HD 26793 (spectral type B9V) and HD 
74280 (B3V), obtained before or after the SN observa- 
tion on the same night with the same setup. The slit of 
the DDO spectrograph has fixed E-W direction, so the 
observations could not be taken at the parallactic angle. 
During the observations the airmass of the SN was always 
less than 2, but differential light loss may have reduced 
the observed fluxes in the blue part of the spectra (espe- 
cially below 4500 A). However, this effect is probably not 
strong, because the shapes of our spectra are similar to 
those collected from literature (see below). 
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Fig. 4. The bolometric light curve of SN 2002ap. The con- 
tinuous Une represents the fitted analytic function given 
by Eq.l, while the dotted line shows the rate of the energy 
input by the indicated radioactive decay. 



Table 3. Journal of spectroscopic observations. The 
coloumns contain the following data: date and JD of the 
observation, phase relative to B-maximum (JD 2452311.5, 
see text), central wavelength (in A), spectral resolution, 
airmass of the SN and the parallactic angle at the time of 
the observation (in degrees). 



Date 


JD-2452000 


Phase 


Ao 


A/AA 


2002.02.02 


308.49 


-3 


5800 


800 


2002.02.03 


309.50 


-2 


5800 


800 


2002.02.06 


312.54 


+1 


5890 


14700 


2002.02.09 


315.54 


+4 


5800 


800 


2002.02.15 


320.53 


-F9 


5800 


800 


2002.02.18 


323.50 


+ 12 


5800 


800 


2002.02.24 


330.49 


+19 


5800 


800 


2002.02.25 


331.52 


+20 


5800 


800 



1.22 
1.25 
1.46 
1.60 
1.54 
1.45 
1.60 
1.82 



After the calibrations, the spectra were de-reddened 
using E{B — V) = 0.09 mag (discussed below) and cor- 
rected for the redshift of the host galaxy (z = 0.0021, 
[Foley et al. 2003] ). Table 3 contains the journal of spec- 
troscopic observations. 

The final, redshift-corrected and dereddened spectra 
are plotted in Fig. 5, together with the spectra taken by 
IGal-Yam et al. 20021 at Wise Observatory. The character- 
istic features of hypernovae are apparent: very blue, al- 
most featureless spectrum at early phases, very broad, 
strong troughs around maximum light, and flat continuum 
with shallow, broad absorption lines at the end of the pho- 
tospheric phase. These spectra are very similar to those 
presented by [Kinugasa et al. 2002| and [Foley et al. 2003| 
Kinugasa et al. showed that the deep absorption troughs 
developing around maximum are due to Can, Fell, Sill, 
Nill, Coll, Ol and Nal. 
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Fig. 6. The velocity curve of hypernovae 
SN 2002ap (filled c ircles), SN 1998bw (open cir- 
cles, IPatat et al. 2001|l and SN 1997ef (open squares, 
[Mazzali, Iwamoto fc Nomoto, 2006) 1. 



Radial velocities were determined from the usual 
Sin A6355 line. Table 4 lists the derived velocities, 
together with other available data collected from the 
literature. The uncertainties were estimated from the 
width of the Sin absorption trough. These errors are 
rnthe r large due to the strong broadening at the very 
early phases and the decreasing line strength at late 



lases. Fig. 6 contains the comparison of Sin veloci- 
t||s of hypernovae SN 1998bw l|Patat et al. 2001|l and 
1997ef dMazzah, Iwamoto fc Nomoto, 2000| ). This fig- 
agrees very well with a similar plot presented by 
[Ijiiugasa et al. 2002[ and supports their conclusion that 
SJS 2002ap evolved faster than the other two well-observed 
J^jernovae (see also [Foley et al. 2003| ). The radial veloc- 



ities will be used in Sect. 3 and 4. 

A single high-resolution spectrum was obtained on 6th 
Feb. 2002 centered around Na D (Fig. 7). The presence of 
the narrow dublet at zero velocity is clearly indicative of 
interstellar absorption originating from Milky Way ISM. 
The equivalent width of the two components was mea- 
su red to be 0.30 ± 0.0 4 and 0.31 ± 0.05 A. According 
to [Barbon et al., 1990[ the total EW of the dublet cor- 
responds to E{B -V) = 0.15 ± 0.02 m ag. The more 
recent relation of [Munari fc Zwitter, 1997] (that uses only 
the EW of the A5891 A component) results in 0.12 ± 0.01 
mag. The uncertainties given here are formal errors of 
the EW measurement errors. A more realistic error bar 
that refiects the uncertainties of the applied empirical 
relations is in the order of ±0.1 mag. Despite of the 
lower S/N of our high- resolution spectrum, these estimates 
agree quite well with the value derived by [Foley et al. 2063| 
{E{B — V)~ 0.1). This confirms the conclusion of many 
previous papers that the reddening of SN 2002ap is low, 
but not negligible. 
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Fig. 5. The spectral evolution of SN 2002ap around maximum (solid line: DDO spectra, dotted line: Wise Obs. 
spectra). The phase relative to _B-maximum is indicated on the right-hand side. Each spectrum was shifted vertically 
for better visibility. 



For such low reddening, the all-sky map of 
ISchlegel et al. 1998| usually works better. This gives 
E{B — V)gai = 0.071 mag as the galactic component 
(that is also probed by the observed zero velocity NaD 
dublet). This is somewhat lower than the values listed 
above, but it is within the global uncertainty of the above 
estimates. The contribution of the host galaxy ISM was 
measured bv lTakada-Hidai et al. 20021 and found to be as 
low as 0.02 mag. The sum of these two components gives 
E{B — V)= 0.09 mag, which was used in almost all pre- 
vious papers on SN 2002ap. Therefore, we have adopted 
E{B -V) = 0.09 ± 0.02 mag (estimated uncertainty) 
throughout the rest of this work. Our spectroscopic re- 
sult, although with higher uncertainty, supports this low 
reddening value. 



3. Distance estimate 

Distance is one of the most important parameter for SNe, 
because the physical parameters, such as luminosity and 
kinetic energy, can be inferred reasonably only if the dis- 
tance is known. Core-collapse events, such as hypernovae, 
usually exhibit large variation in their peak luminosity and 
ejected mass | |Hamuy, 2003| l, thus, distance determination 
is necessary to infer these parameters from the observed 
ones. 

The distance of a hypernova is not an easy parame- 
ter to determine. Because these events are far from being 
"standard candles" , a pure photometric distance measure- 
ment is not possible at present. A possible solution for this 
problem may be the determination of the distance of the 
host galaxy itself, via another, independent method. This 
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Table 4. Radial velocities of SN 2002ap inferred from the 
Doppler-shift of the Si ii A6355 line. The second column 
gives the phase with respect to the date of explosion (JD 
2452303.4. I Mazzali et al. 2002|l . The errors of the radial 
velocities (fourth coloumn) were estimated from the width 
of the Sin line. The meaning of the reference codes is: 
K - IKinugasa et al. 2002| W - IGal-Yam et al. 20021 pp - 
present paper. 



JD t t;,.(Si II) cr ref. 

(2452000+) (days) (kms~^) (kms^^) 

K" 
W 
K 
PP 
PP 
K 
K 
W 
K 

PP 
K 
K 
K 
K 
W 

PP 
K 

PP 
K 
K 
W 
K 

PP 
K 

PP 



might work at least for the nearby galaxies, such as M74. 
Surprisingly, such pre-determined distances for M74 are 
quite uncertain (see Section 1). Thus, in order to obtain 
reliable parameters, a distance measurement of SN 2002ap 
itself would be very useful. 

3.1. The Expanding Photosphere Method 

We have applied the Expanding Photosphere Method 
(EPM) for measuring the distance to SN 2002ap. This 
method was developed by [Kirshner fc Kwan, 1974| and 
was used many times for measuring distances of Type 
II SNe fe.g- .Schmidt et al. 1994l|Hamuy et al. 2001| ). The 
main assumptions of this method are as follows: 

— the expansion is spherically symmetric and homolo- 
gous 

— there is a well-defined photosphere in the ejecta, below 
which the gas is optically thick 

— the photosphere radiates as a (diluted) blackbody 
(parametrized by the dilution factor Q. 
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Fig. 7. Na D fines in the spectrum of of SN 2002ap. 



Because Type II SNe eject a massive H-rich envelope, 
they have optically thick atmospheres, in which the loca- 
tion of the photosphere is determined by the position of 
the H recombination front. This phase lasts for months, 
therefore the EPM can be applied for data obtained over 
a long time interval. 

It is very probable that none of the above assumptions 
are fully applicable for SN 2002ap, which was very dif- 
ferent from a Type II SN. According to the spectroscopic 
observations (including those presented in this paper), the 
photosphere disappeared quickly, and after a few weeks 
the SN entered the nebular phase. Therefore, the EPM 
can be applied only for a few data points just before and 
after maximum light. Moreover, the atmosphere of hy- 
pernovae could be very different from that of a Type II 
SN, and the emergent radiation may be far from being 
Planckian. In their very thorough study of Type II SN at- 
mospheres, Eastman ct al. 1996 showed that the dilution 
factor C for a particular atmosphere is usually between 
0.4 - 1.5 depending on the effective temperature, density 
profile and the choice of the wavelength regime for deter- 
mining the colour temperature. Because Q directly influ- 
ences the EPM-based distances, neglecting the deviation 
from Planckian flux can lead to a 40% systematic error in 
the derived distance. For hypernovae, such kind of prob- 
lems can be studied only via detailed model atmospheres, 
which is beyond the scope of this paper. Nevertheless, we 
assumed that the EPM can be applied to SN 2002ap, as 
a first approximation, until the more sophisticated model 
atmosphere computations are available. 

Perhaps the weakest point of the feasibility of the 
EPM is the assumption of spherical symmetric expansion. 
According to the polarization results (Section 1), there are 
indications that the expansion of SN 2002ap was definitely 
asymmetric at early phases. The measured net polariza- 
tion (about 2 %) may suggest an asymmetry of about 10 
% (|Kawabata et al. 2002|l . although the polarization an- 
gle varied in time, which means that the asymmetry, if 
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indeed present, may be of complex nature. Another pos- 
sibility may be the presence of higher interstellar polar- 
ization than Kawabata et al. estimated. In this case the 
measured polarization is not due to the elongated shape of 
the ejecta. Indeed, recently |Leonard"et al. 2002b has dis- 
covered evidence for extremely high polarization efficiency 
in NGC 3184, where a high degree of polarization was 
measured for SN 1999gi that otherwise was only slightly 
reddened. It is also worth noting that model computations 
for SN 2002ap could reproduce the early part of the light 
curve using sperical as well as non-spherical models (see 
e.g. lMaeda et al. 2003l for details). In the followings we as- 
sume, for simplicity, that the expansion of SN 2002ap was 
spherically symmetric during the photospheric epoch, but 
note that deviations from sphericity may result in sys- 
tematic errors in the distance determination. Discussion 
of other possible sources of random and systematic uncer- 
tainty is presented at the end of this Section. 

The formulae of the EPM are straightforward. 
According to the homologous expansion, the radius of the 
photosphere can be expressed as: 



R — Vphot{t — to) + i?0 



(2) 



where is the moment of explosion, to is a crucial pa- 
rameter for the distance determination. The angular size 
of the photosphere, 9 = R/D, where D is the distance, 
can be expressed as ( |Kirshner fc Kwan, 1974| ) 



fx 



C>SA(rA) 



(3) 



where f\ is the observed flux density, B\{Tx) is the Planck 
function and (^\ is the dilution factor. 

In this paper we used the UVOIR bolometric light 
curve of SN 2002ap (Section 2) for estimating the surface 
brightness - temperature relation. The dilution factor is 
assumed to be C = 1 (e^s a first approximation) and Eq.3 
can be written as 



/ fbol 



(4) 



where fboi is the observed bolometric flux density (in 
erg/s/cm^) and a is the Stefan-Boltzmann constant, fboi is 
given from the bolometric light curve discussed in Sect. 2. 
The effective temperature was estimated from hypernova 
model atmospheres. 

We used the hypernova model atmospheres by 
[Mazzali, Iwamoto &: Nomoto, 2000| constructed for 
SN 1997ef. For several epochs, they tabulated the effec- 
tive temperature and the calculated B — V color of the 
photosphere. We fitted a simple linear function to these 
data and found 



T, 



0.122(±0.03) (S - y) -I- 3.875(±0.03). (5) 



Combining Eq.4 and 5, the angular size of the photosphere 
9 can be derived from the observations. 
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Fig. 8. The effective temperature of SN 2002ap (via Eq.5) 
as a function of time. The filled circles correspond to data 
obtained during the expanding photosphere phase (see 
Fig.9). 
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Fig. 9. The variation of the angular size (in mas) of the 
photosphere of SN 2002ap. The expansion of the photo- 
sphere lasted only about 12-13 days after explosion. 



The photospheric velocities were assumed to 
be equal to the radial velocities based on the 
usual Sin line (Sect. 2). The hypernova models of 
[Mazzali, Iwamoto fc Nomoto, 2000] give some support 
to our assumption, as their computed Si ii velocities are 
similar to the photospheric ones. 

Having the angular size and the velocity determined 
from the observations at several epochs, the distance can 
be derived via Eq.2 and 4. Following .Schmidt et al. 19941 
their combination yields 



t ^ D k ( ) + to, 

Vphot 



(6) 



where D is in Mpc, Vphot is in kms^^, t is in days and 
k — 3.5697 • 10^'' is the conversion factor between the 
quantities. Plotting t versus 9/vphot, D and to can be de- 
termined by linear regression. 
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Fig. 10. The fit of Eq.6 to the data of SN 2002ap, yielding 
the distance and the moment of explosion. 



The results for SN 2002ap are plotted in Figs.8, 9 and 
10. Fig. 8 contains the plot of the temperature variation 
(from dereddened [B — V)o data, via Eq.5). Fig. 9 shows 
the angular size calculated from Eq.4. It is apparent that 
the expansion of the photosphere stopped at about JD 
2452317, so the expansion phase did not last long. Note, 
that, as the referee pointed out, this does not mean that 
the photospheric phase has ended, because the spectra 
clearly indicate that the SN is still in the photospheric 
phase (see Fig. 5). The EPM was applied only for data ob- 
tained before this epoch. The fitting of Eq.6 to the com- 
bined data is presented in Fig. 10. The linear regression 
resulted in D = 6.7 ±0.5 Mpc and to = 2452304.3 ±1.0 as 
the distance and the moment of explosion, respectively. It 
is clear from Fig. 10 that the relation becomes non- linear 
when the expansion phase ends, thus, at least one of the 
assumptions of the EPM became invalid at that time. The 
quoted uncertainties are formal (random) errors due to the 
fitting procedure. 

These results show remarkably good agreement with 
previous ones, determined via other methods. The EPM- 
distance is almost perfectly equal to the photometric dis- 
tance of M74 when the |Schlegel et al. 1998| reddening cor- 
rection is apphed {d = 6.8 Mpc, see Sect.l). The moment 
of explosion agrees within the uncertainties with the one 
derived bv IMazzah et al. 20021 (tn = 2452303.4, Jan. 28, 
2002) from model computations. It also worths noting that 
the derived explosion date is only 1 day prior to the dis- 
covery of SN 2002ap (see Sect.l). 

3.2. Uncertainties of the EPM-distance 

As mentioned above, difference between the assumptions 
of the EPM and the true physical state of the SN at- 
mosphere may lead to significant uncertainty in the in- 
ferred distance. Despite the good agreement between 
the new EPM-distance and other distance estimates, 
all these may suffer from quite large systematic errors. 



Recently, ILeonard et al. 20031 presented a thorough test 
of the EPM- and Cepheid-based distances of SN 1999em 
and NGC 1637. They concluded that the Cepheid-based 
distance is 50 % greater than the EPM-based distance, 
which means that there is serious disagreement between 
these two measurement methods. They also showed that, 
at present, the Cepheid-based distance is more favourable 
than the EPM-based one, thus, the EPM-distances are 
plagued by at least 50% systematic error. This leads to a 
60 - 70 % total (systematic plus random) uncertainty for 
the EPM. 

What are the possible sources of such high uncer- 
tainty? ILeonard et al. 20031 mention the dilution factor 
and the radial velocity measurement technique as the 
most probable ones. The need for a precise dilution factor 
based on model atmosphere tuned for the particular SN 
has been emphasized by Eastman et al. 19961 However, 
their dilution factors may also suffer from uncertainies, be- 
cause independent calculations by other groups have lead 
to significantly different factors (see ^Leonard et al. 20031 
for extensive discussion and references). These problems 
are probably more enhanced for the present study, due 
to the atmospheric properties of hypernovae. Since there 
was no detailed model atmosphere for SN 2002ap at our 
disposal, we had to assume C = 1, which may lead up 
to 50 % systematic error in the distance (as it was men- 
tioned in Section 3.1). This error can be even larger if 
the atmosphere of SN 2002ap were very different from 
a Type II SN atmosphere so that the dilution factor 
was quite different from those of Eastman et al. 19961 
On the other hand, the application of a hypernova 
model of Mazzali, Iwamoto fc Nomoto, 2000| (constructed 
for SN 1997ef) may somewhat reduce this large uncer- 
tainty. As it is explained above, we assigned the model 
effective temperature (via Eq.5, by comparing the ob- 
served and predicted B — V color indices) to the bolometric 
flux data. This technique may have the advantage of us- 
ing the UVOIR "bolometric" fluxes that are integrated 
over a long wavelength regime, thus, it is not sensitive 
to the chosen wavelength and wavelength-dependent dilu- 
tion factors. On the other hand, the uncertainties of the 
derived bolometric fluxes are probably in the same order 
of magnitude as that of the dilution factor. 

The systematic errors caused by the measurement 
method of radial velocities are also quite large for 
SN 2002ap because of the very broad lines and strong 
blends in the spectra, especially at the early phases. 
Therefore, referring to ILeonard et al. 20031 we estimate 
the systematic uncertainty of the EPM-based distance as 
60 % that means ±4 Mpc for SN 2002ap. The total (ran- 
dom plus systematic) error of the distance becomes ±4.5 
Mpc, or about 70 %. 

The consistency between the EPM-distance and 
the other photometric distances of M74 may sup- 
port the applicability of the EPM for hypernovae. 
However, it is important to note that ILeonard et a l. 2003j 
found similar agreement between the EPM-distance 
and the "brightest red supergiants" (BRSG) dis- 
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tance dSohn fc D avidge, 1998| ) for NGC 1637. Among 
the distance measurement methods discussed by 
ILeonard et al. 20031 the BRSG technique is the only 
one that gives consistent result with the EPM. Therefore, 
this agreement does not, unfortunately, reduce the large 
systematic uncertainty of the EPM-distance for M74. 

We conclude that our first attempt to apply the EPM- 
technique to a hypernova resulted in distance estimate for 
SN 2002ap that has remarkably good internal precision 
(good linearity, low scattering in Fig. 10.) On the other 
hand, the systematic uncertainty of the EPM-distance is 
quite high due to the problematic issue of the dilution fac- 
tors as well as the complexity of hypernova spectra at the 
early phases. Note that the EPM-distance is 0.6 Mpc lower 
than the one used in previous papers (7.3 Mpc), thus, 
the inferred physical parameters of SN 2002ap may need 
some revision, if we adopt this new distance. For exam- 
ple, the nickel mass may be somewhat lower than derived 
previously (e.g. IMazzah et al. 20021 [Pandey et al. 2002| 
[Foley et al. 2003| ). 



4. Model computations and physical parameters 

Some basic physical parameters of the SN explosions can 
be estimated by comparing the observations with model 
computations of the light curves and/or spectra. This usu- 
ally requires sophisticated calculations including hydrody- 
namics, radiative transfer and atomic physics, which are 
beyond the scope of this paper. Instead, in the foUowings 
we use a simple analytic ("toy"-) model to estimate the 
fundamental parameters of SN 2002ap, such as Ni-mass 
and kinetic energy of the ejecta. 

The bolometric light variation of SNe can be qualita- 
tively described by simple analytic m odels l| Arnett, 1980| 
jArnett, 1982| llwamoto et al. 20021 IMaeda et al. 2003|l ' 
assuming homologous expansion of the ejected envelope 
and (probably over-) simplified deposition of gamma-rays. 
It is known that SNe are powered by radioactive decay 
of ^^Ni — >• ^^Co ^^Fe. The basic parameter characteris- 
ing the emitted energy is Mat^, the mass of ^^Ni synthe- 
sized during the explosion. The gamma-rays originating 
from radioactive decay are deposited and thermalized in 
the ejecta via Compton-scattering. Because the ejecta is 
optically thick at early phases, the timescale of the radia- 
tive transfer of photons to the surface is comparable to 
the timescale of the expansion (characterized by v^xpi the 
velocity of the outmost part of the ejecta), thus, the en- 
ergy remains trapped in the optically thick atmosphere. 
As the ejecta expands, the density and the opacity de- 
creases, and the deposited energy can escape more and 
more quickly. Thus, the light curve has a peak at 15-20 
days past explosion, when the emitted luminosity becomes 
roughly equal to the instantaneous rate of the energy de- 
position ( |Arnett, 1982| ) . Because the transparency of the 
ejecta quickly decreases with time, the light variation af- 
ter maximum mimics the radioctive decay law and the 
physics of gamma-ray deposition. 



The very simple light curve model used in this paper 
is based on the assumption of free, homologous expansion 
of the ejecta. The ejecta has a radius R{t) that increases 
with time as 



R{t) = i?o 



(7) 



where i?o is the (negligible) initial radius of the progenitor 
and Vexp is the (constant) expansion velocity. The velocity 
of a thin shell at fractional radius x = r{t)/R{t) has a 
velocity 



v{x) 



^exp • 



(8) 



The density structure of the ejecta is assumed to be 
consisting of an inner core with constant density, and an 
outer shell where the density is decreasing as power-law: 



(9) 



where xq is the fractional radius of the core, po is the 
density in the core (i.e. for x < xg) and n is the power-law 
exponent. As a consequence of the expansion, po decreases 
in time as 



Poit) 



M 



1 



(10) 



where /(xq) = (xq^/S -|- a;o"(l — xo^~")/(3 — n)) is a ge- 
ometric factor due to the assumed density structure, and 
M is the total mass of the ejecta. With these assumptions 
the kinetic energy can be expressed as 



Eh — — 



3 3 — n 5xn 



10 5-n 3a;" 



5 

—I ^^-p 



nx: 



(11) 







The energy production released as gamma-rays is de- 
scribed by the exponental decay law: 



= Mm {Sm + 0.82 Sco) 
Sni — Effi exp(— Aatj^) 
Sco = Eco (exp(-Acoi) - exp(-AArit)) 



(12) 
(13) 
(14) 



where Eni = 3.70 • 10^° erg/g, Eco = 6.76 • 10^ erg/g, 



Ajvi = 1.3152 • 10-^ s-i, Aco = 1.0325 • 10 



-7 



, and 



M^i is the total mass of the synthesized ^^Ni. A fraction 
of the decay energy goes to generating positrons, which 
are decelerated and annihilated by the ejecta (see e.g. 
jCapellaro et al. 1997| ). The energy input due to positrons 
is approximated as 



E^ 



0.08 MN^ Sco 



(15) 



Thus, the total energy generation rate becomes the sum 
of the gamma ray and positron energy productions. 

The gamma-rays as well as positrons must interact 
with the ejecta through multiple scattering in order to 
deposit their energy, which is assumed to be fully thermal- 
ized and re-radiated as optical photons. The total amount 
of deposited energy at each time t is approximated as 



e = E^ (1 - exp(-r^)) + E+ {1 ~ exp(-r+)) 



(16) 
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Table 5. Computed models of SN 2002ap. The meaning of the symbols are: M: ejected mass (in Mq); M^i'- mass 
of synthesized Ni (in Mq); xq: fractional core radius; n: power-law exponent of density; Vexp'- expansion velocity (in 
(kms^^); K-yi 7-ray opacity (in cm^g~^); Kg+i positron opacity; Kopt- opacity for optical photons; -Esi: kinetic energy 
of the explosion (in 10^^ ergs) The asterisk indicates that the parameter was kept fixed during the fitting. 



Model 


M{Mq) 


A'/iv^M©) 




n 








f^opt 


-E51 


A 


1.0 


0.07 


1.0* 


0.0* 


35000 


0.1 


5.0 


0.009 


7.35 


B 


1.0 


0.04 


0.01* 


2.1 


35000 


0.005 


0.1 


0.03 


3.84 


CI 


1.0 


0.07 


0.15* 


2.1 


40000 


0.027* 


1.4 


0.025 


5.67 


C2 


2.0* 


0.045 


0.15* 


2.1 


35000 


0.027* 


1.0 


0.012 


8.69 



where 



1 — X, 



n-l 



= Po(0 R{t) Xq (IH ) 

n — 1 



(17) 



is the total optical depth for gamma-rays and positrons 
( |Capellaro et ai. 1997| ). The opacity k, is assumed to be 
constant for both gamma-rays and positrons. This is a 
major simplification, but often used in SNe models (e.g. 
ICapellaro et al. 1997| Pinto fc Eastman, 2000D . 

If the diffusion time is short, the emitted energy per 
second is approximately equal to the rate of the de- 
posited energy, i.e. Lopt ~ e. We have adopted this 
approximation for calculating the bolometric luminosity 
(similarly to IMaeda et al. 200311 . instead of solving the 
time-dependent transport problem for optical photons 
(see |Capellaro et al. 1997| |. This condition is certainly not 
valid for pre-maximum epochs, thus, our simple light curve 
model is applicable only on the declining part of the light 
curve. 

The light curve model outlined above contains the 
following free parameters: the total ejected mass M, 
the nickel mass Mjvi, the expansion velocity v^xp, the 
density exponent n, and the opacities K-y and k+. 
These parameters can be constrained by the observa- 
tions. However, it is known that the light curve alone 
does not determine all parameters uniquely, i.e. the same 
light curve could be described by several parameter 
combinations Ijiwamoto et al. 19981 llwamoto et al. 2002|l . 
Therefore, additional (usually spectroscopic) data are 
needed to constrain the physical parameters. Since our 
code is not capable of modeling the spectrum, we have 
used the Si 11 velocities for this purpose. 

In a grey atmosphere (such as our simple model with 
constant opacity) the location of the photosphere Xph is 
usually defined as 



Kopt p{.x,t) R{t) dx 



(18) 



Using Eq.2, the photospheric velocity can be expressed as 



Vph 



^exp 



[1 



2(1 -n) 



1 



ll/l-n 



«opt Poit) x'S R{t) 



(19) 



We have fitted Eq.l6 and 19 to the bolometric light 
curve and the radial velocity curve presented in Sect. 2. 



Three kinds of models were defined depending on the den- 
sity configuration (characterized by the density param- 
eters Xq and n): a constant-density model with n — 
fixed (Model A), a power-law model with xq = 0.01 fixed 
(Model B) and a "core-shell" model with xq = 0.15 fixed 
(Model C). In the latter two, the density exponent n was 
treated as variable. In addition, the parameters M, M^i, 
Vexp, ^^7, and Kopt were varied to find the best agree- 
ment between the observed and computed data. 

The fitted parameters are collected in Table 5. In gen- 
eral, low-mass models with M = 1 Mq gave adequate fit 
to the observed data. All models have shallow density de- 
pendence (n w 2). The kinetic energy of the explosion E^i 
is between 4-9 and the nickel mass is Mjvi = 0.04 - 0.07 
Mq. For Model C, we have also forced an M — 2 Mq 
model fit (model C2 in Table 5). This model resulted in 
slightly worse, but still acceptable fit (with less Ni mass, 
expansion velocity and opacity values). Similar tendency 
was found in the case of Models A and B: the increase 
of the ejected mass required the decrease of the Ni mass, 
velocity and opacities. It indicates that the parameter de- 
generacy, mentioned above, is not completely solved by 
introducing the velocity data. 

Fig. 11 and 12 show the comparison of Model CI (solid 
line) and C2 (dashed line) with the observations. Model 
A and B resulted in very similar light and velocity curves 
(note that the models are valid only for the late-time light 
curve) . 

Because only Model C could be fitted with the usual 
value of the gamma-opacity = 0.027 cm^g"^, we have 
adopted model CI as the most probable. This model is 
somewhat similar to the two-component model described 
by IMaeda et al. 200'3l because the density structure and 
the optical depth is different for the core and the outer 
shell. IMaeda et al. 20081 discussed that a dense core is re- 
sponsible for the late-time decline of the light curve. Note 
that our simple constant density model (Model A) can 
also describe the late-time light curve, but this requires a 
radical increase of the gamma-ray and positron opacities, 
which is difficult to explain physically. 

Indications for a shallow density profile of the 
ejecta was also given by spectral synthesis models. 
Using the parametrized spectral synthesis code SYNOW, 
IKinugasa et al. 2002| obtained r ~ {v/vph)^^ for the op- 
tical depth above the photosphere. Assuming homologous 
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Table 6. Comparison of model parameters of SN 2002ap. 40000 
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Fig. 11. Comparison of the computed bolometric light 
curves with the observations. Continuous line: Model CI; 
dashed line: Model C2 (see Table 5 for the list of param- 
eters). 



expansion (v ^ x) 
corresponds to p 



ifMazzali ct al. 2002|l also resulted in a flat density distri- 
bution above the photosphere. The result of our simple 
model {p ~ x~'^) is consistent with these more sophisti- 
cated calculations. 

Table 6 summarizes the 



{Mej,Mm,E5i) of SN 2002ap 



basic parameters 
, jv^ivi, j^5i; <~>^i ^uu^ay coUectcd from the 
literature. It is seen that our Model CI has lower ejected 
mass than the previous models, but the Ni-mass and the 
kinetic energy is similar. This is a bit surprising, because 
we used a somewhat shorter distance modulus than the 
authors in the previous papers, thus, our bolometric light 
curve is slightly fainter. Because the peak luminosity is 
directly connected with the Ni-mass, one should expect 
a lower Ni-mass for a fainter peak luminosity. However, 
our light curve code is probably too simplistic for such 
fine-tuning of the physical parameters of SN 2002ap, thus, 
the parameters listed in Table 6 can be considered only 
as order-of-magnitude estimates. As IMaeda et al. 20031 
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Fig. 12. Observed (dots) and computed photospheric ra- 
dial velocities. The line types correspond to the same mod- 
els as in Fig. 11. 



pointed out, a two-component (or core-shell) ejecta model 
may, in fact, be an indication of an asymmetric explosion. 
If this is indeed the case, the parameters derived from 
such spherical models may deviate significantly from the 
real values. 

We conclude that the simple model of the light- and 
velocity curves of SN 2002ap confirms that the progenitor 
was a low- mass object (probably a stripped stellar core), 
and the synthesized nickel mass was quite low, about 0.07 
Mq, or even less, depending on the density structure of 
the ejected envelope. 

5. Summary 

The results of this paper can be summarized as follows. 

1. We have presented new BVRI photometry of SN 
2002ap up to t = 300 days past maximum, and photo- 
spheric phase optical spectroscopy. 
The Expanding Photosphere Method was used to in- 
fer the distance of SN 2002ap. This resulted in d = 
6.7 Mpc (corresponding to /io — 29.13 mag) with 
±0.5 Mpc random and ±4 Mpc (estimated) system- 
atic uncertainty. This agrees with previous photo- 
metric distances of M74, if the reddening map of 
[Schlegel et al. 1998| is used to estimate the galactic 
reddening. However, it is slightly less than the value 
(7.3 Mpc) used in previous papers of SN 2002ap. 
The physical parameters of the SN ejecta were esti- 
mated with a simple parametrized analytic model, as- 
suming homologous expansion and power-law density 
structure. From simultaneous fitting to the bolometric 
light curve and the radial velocity curve, we found that 
the SN ejected a low mass envelope {M^j « I.OMq) 
and synthesized about 0.04 - 0.07 Mq of ^^Ni. 
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